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ABSTRACT 


Adaptive  control  system  techniques  have  been  utilized  to  inves 
tigate  identification  of  the  dynamic  equations  of  operating  control 
systems.  The  identification  has  been  achieved  by  using  model  modu 
lation  techniques  to  drive  dynamic  models  into  correspondence  with 
operating  control  systems.  The  system  identification  then  proceeded 
from  examination  of  the  model  and  the  adaptive  loop. 

The  model  modulation  techniques  applied  to  adaptive  control 
systems  are  breifly  discussed  These  techniques  applied  to  simulation 
studies  on  an  analog  computer  and  a  digital  computer,  and  applied 
directly  to  a  position  servo-mechanism ,  are  then  discussed  The 
discussion  presents  data  showing  the  accuracy  of  identification  of 
the  dominant  characteristics  of  the  dynamic  equations  when  the  order 
of  the  equation  is  either  known  or  approximated  It  was  found  that 
these  dominant  roots  could  be  found  within  an  average  accuracy  of  570 
for  both  complex  phase  angle  and  undamped  natural  frequency 

The  authors  wish  to  express  their  appreciation  to  Professor 
Richard  Co  Dorf  of  the  U.  So  Naval  Postgraduate  School,  for  the 
suggestion  of  the  model  modulation  method  and  for  his  assistance 
during  this  investigation. 
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CONTROL  SYSTEM  IDENTIFICATION  THROUGH  MODEL  MODULATION  METHODS 


Chapter  I 
INTRODUCTION 


Genera  1 

The  last  three  decades  have  seen  control  system  technology  grow 
from  a  series  of  rules  of  application  to  a  full  regime  encompassing 
most  of  the  disciplines  of  engineering.  The  study  of  control  systems 
has  progressed  deeply  into  analysis  and  synthesis  of  linear  and  non¬ 
linear  control  systems.  There  are,  additionally,  offshoots  of  basic 
systems  with  most  exhaustive  methodology:  examples  are  the  sampled 
data,  hybrid  analog-digital  and  multi-loop  self  adaptive  systems. 

The  problem  of  control  system  identification  led  to  some  of 
the  original  investigation  in  the  control  field.  Despite  the  increase 
of  knowledge  in  the  field  of  control  engineering,  system  identification 
has  remained  a  major  problem. 

The  intent  of  this  project  was  to  investigate  a  method  of  control 
system  identification  utilizing  techniques  from  the  growing  technology 
of  the  self-adaptive  control  system.  One  restriction  was  laid  upon 
the  project  at  its  inception.  This  was  that  the  method  of  system 
identification  remain  independent  of  the  system  to  be  identified 
Ideally  then,  an  operating  system  should  be  identified  solely  through 
examination  of  the  input  and  the  response.  To  this  end  a  control 
system  identifier  has  been  modeled  and  tested  on  an  analog  computer, 
through  digital  computer  analysis,  and  upon  a  position  control  system. 


Control  System  Identification 

The  characteristics  and  performances  of  the  fundamental  components 
of  feedback  control  systems  are  generally  amenable  to  systematic  analysis 
The  control  system,  when  composed  of  these  fundamental  components,  may 
be  easily  synthesized  and  analyzed.  Components,  however,  when  removed 
from  the  theoretical  state  do  not  exhibit  the  invariant  characteristics 
attributed  to  them.  Furthermore,  control  systems  are  often  synthesized 
on  the  basis  of  transfer  function  manipulation.  This  synthesis  occurs 
under  the  assumption  of  negligible  loading  of  one  component  system  by 
another.  This  assumption  is  too  often  fallacious.  In  either  or  both 
of  these  situations  the  control  system  engineer  is  confronted  with  a 
control  system  which  does  not  conform  to  previous  analysis.  And  if 
this  analyzed  control  system  is  a  portion  of  some  larger  process,  it 
may  be  necessary  to  reevaluate  the  control  characteristics  and  perform¬ 
ance,  in  other  words  to  identify  the  system. 


In  some  extremely  variant  systems  the  processes  of  system  identi¬ 
fication  are  utilized  to  serve  as  controllers  in  a  self-adaptive  control 
system.  The  most  notable  of  these  cases  are  the  self-adaptive  control 
systems  currently  in  use  and  being  further  developed  for  aircraft 
installation.  The  NASA  X-15  is  an  example.  In  this  installation  the 
response  characteristics  are  examined  and  the  results  are  used  to  vary 
flight  control  loop  gain  to  maintain  near  invariant  control  response. 

The  equipment  used  for  this  process  is  quite  complex,  but  is  justified 
by  the  operating  environment  of  the  vehicle. 
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For  less  complex  control  systems  or  in  less  complex  environments, 
the  problem  of  control  can  be  solved  by  more  conventional  means,  without 
recourse  to  dynamic  self-adaptive  system  technology.  These  less  complex 
systems  require  accurate  identification  of  the  process  dynamics.  Once 
the  transfer  function  of  a  process  has  been  evaluated  this  evaluation 
can  serve  as  a  basis  for  design  of  the  remainder  of  a  system  or  for 
readjustment  of  the  process  through  compensation  techniques. 

Identification  Methods 

Truxal  writing  in  reference  (1)  has  placed  the  methods  of  control 
system  identification  into  four  broad  categories: 

a.  Parameter  identification. 

b.  Frequency  response  methods. 

c.  Impulse  response  methods. 

d.  Direct  evaluation  of  system  differential  equations 
of  motion. 

Some  of  the  difficulties  in  one-by-one  determination  of  model 
parameters  have  been  mentioned.  Additionally,  this  method  is  possible 
only  when  the  detailed  operation  of  a  process  is  reducible  to  funda¬ 
mental  differential  equations  in  which  the  primary  parameters  may  be 
included  and  the  negligible  ones  discarded.  This  ideal  situation 
rarely  exists  for  complex  control  of  complex  processes. 

As  frequency  response  techniques  are  well-known,  a  description 
of  them  has  been  omitted.  These  procedures  are  extremely  important 
in  general  control  design.  This  is  the  most  common  method  of  iden¬ 
tifying  the  transmittive  dynamics  of  a  control  process.  Many  working 
control  systems  however,  cannot  be  subjected  to  variable  frequency 
oscillations  for  reasons  of  system  complexity,  operating  economy,  or 
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because  of  the  extensive  time  and  effort  required  for  determination 
of  the  gain  and  phase  characteristics..  Two  additional  problems  occur. 
First,  in  some  circumstances,  frequency  response  testing  should  not 
be  superimposed  on  regular  operation,  because  testing  would  interfere 
with  designed  operation.  Secondly,  a  working  system  is  subjected  not 
only  to  its  driving  input,  but  also  to  noise  occuring  during  the  op¬ 
erating  process.  The  extraneous  disturbances  of  normal  operation 
would  make  the  interpretation  of  frequency  test  results  difficult. 


The  impulse  response  can  serve  to  completely  identify  any  linear 
time-invariant  two  port  system.  The  various  impulse  response  techniques 
are  utilized  in  many  of  the  more  sophisticated  adaptive  control  systems. 
Truxal  has  cited  three  primary  techniques  for  measuring  impulse  response 
The  first  and  most  direct  method  is  to  excite  a  system  with  an  impulse 
and  measure  the  response.  Second,  the  system  can  be  excited  by  a 
known  signal.  For  an  initially  inert  system,  the  impulse  response 
g^t^  can  be  determined  from  the  input,  r^^,  and  output,  g^^,  through 
the  convolution  integral 

A 

I  11 

Finally  the  impulse  response  can  be  measured  by  application  of  a  small 
random  signal  to  the  process.  If  white  noise  is  used  as  the  input, 
the  cross  correlation  function  between  the  output  and  the  noise  is  the 
impulse  response. 


c(t)  I  J(r) r  (t  -y)  ^ 


The  attempt  to  utilize  the  impulse  response  in  identification  is 
complicated  by  three  problems.  Noise  has  a  much  more  delerterious 
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effect  on  the  output  than  in  the  case  of  frequency  response  tests.  Non- 
linearities  are  much  more  difficult  to  detect  than  in  frequency  response 
tests.  Finally  the  mechanization  necessary  to  determine  the  impulse 
response  is  almost  prohibitive  for  use  in  system  identification.  When 
the  impulse  response  method  is  used  in  adaptive  system  control  the 
engineer  is  generally  interested  not  in  the  full  identification,  but 
in  identifying  only  some  characteristic  of  the  response.  Quantitative 
measurements  are  made  to  serve  as  sufficient  identification  to  allow 
adaptive  correction.  As  an  example,  the  number  of  oscillations  during 
a  specific  period  of  a  system  response  would  identify  the  impulse 
response.  These  approximations  are  not  sufficient  for  full  identifi¬ 
cation,  however. 


The  evaluation  of  the  differential  equation  for  a  linear  process 
offers  an  additional  method  of  identification.  The  differential  equation 
of  a  process  might  be  expressed  as  an  equation  of  input  r^^  and  output 
c^^  in  the  following  form: 


& +  a~jf*  -  *  +  +  •••*  b‘r° 


The  a^  term  may  be  evaluated  from  the  nature  of  the  zero  frequency  or 
steady  state  gain.  If  the  order  of  the  differential  equation  is  known, 
and  the  nature  of  the  right  hand  side  is  known,  the  other  ^  co¬ 
efficients  may  be  evaluated  by  generalizations  such  as  minimizations 
of  square  of  the  error  e,  where  e  is  defined: 


e 


TtN 


4 


o  o  0  4 


c,'jr  * 


q  c 


1.3 
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The  ak,  -  coefficients  are  evaluated  through  use  of  this  error  criterion 

N-  1 

by  generation  of  successive  differentiations  of  the  response.  Unfor¬ 
tunately,  searching  over  an  (N-l)  space  for  a  minimum  is  a  problem  of 
potential  difficulty.  General  problems  also  arise  in  settling  at 
local  minimums  and  in  error  system  stability. 

Many  of  the  problems  of  systematic  determination  of  the  differ¬ 
ential  equation  coefficients  are  alleviated  by  the  use  of  system-model 
reference.  This  is  a  natural  variant  to  the  adaptive  control  system 
technology.  The  exploitation  of  this  method  involves  utilization  of 
the  difference  between  system  and  model  response  to  the  forcing  input. 

This  difference  has  been  handled  by  several  methods  in  present  technology. 
In  general  self  adaptive  processes,  the  difference  as  operated  upon 
by  an  adaptive  technique  is  used  to  drive  the  system  into  operating 
correspondence  with  the  invariant  response  of  the  model. 

This  investigation  was  based  on  a  variation  of  this  method.  During 
this  investigation  adaptive  techniques  have  been  used  to  drive  the 
model  into  operating  correspondence  with  the  operating  system.  The 
characteristics  of  the  model  differential  equation  have  been  modified 
by  the  adaptive  loop.  These  characteristics  then  have  been  used  to 
accurately  define  the  coefficients  of  the  system  differential  equation. 
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Chapter  II 


MODEL  REFERENCE  ADAPTIVE  SYSTEMS 


Introduction 

A  major  advantage  of  the  model  reference  adaptive  system  is  that 
direct  measurement  of  the  operating  system  differential  equation 
parameters  is  not  required.  This  factor  results  in  considerable 
economy  in  mechanization  of  the  adaptive  components.  As  previously 
mentioned  there  are  in  existence  various  methods  of  self  adaptive 
control.  However,  in  many  of  these  systems  the  adaptive  loop  mech¬ 
anization  requires  the  assumption  of  the  second  order  nature  of  the 
system.  An  example  of  this  is  Osder's  system  detailed  in  reference  (1). 
For  this  system,  the  sign  reversals  in  the  output  oscillations  are 
counted  and  the  count  used  to  correct  the  second  order  zeta.  Most  of 
the  impulse  response  self  adaptive  systems  also  require  assumption  of 
a  basically  second  order  system  to  keep  the  equipment  within  allowable 
limits.  The  model  reference  system  relieves  this  restriction.  The 
model  reference  adaptive  system  can  be  utilized  for  higher  order  systems 
about  which  the  second  order  assumption  is  invalid.  Additionally, 
current  literature  presents  several  examples  in  which  model  reference 
allows  the  control  of  several  parameters  within  the  operating  system. 

References  (1)  and  (2)  present  several  methods  of  model  reference 
adaptive  control  systems  which  have  been  proposed  or  developed.  The 
variants  of  these  systems  include  the  Ham-Lang  conditional  control 
system,  which  in  analysis  is  passive  adaptation  through  model  reference. 
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The  other  extreme  is  perhaps,  Braun's  method,  which  utilizes  an  ortho¬ 
normal  spectrum  analyzer  to  solve  the  identification  problem. 

Parameter  Perturbation 

Perhaps  the  most  straightforward  adaptive  mechanization  results 
from  the  use  of  model  reference  with  perturbation  techniques.  McGrath, 
Rajaraman  and  Rideout  in  reference  (3)  have  proposed  and  examined 
various  parameter  perturbation  model  reference  schemes.  Figure  1, 
taken  from  their  discussion  presents  two  variants  of  the  parameter 
perturbation  scheme . 


Fig.  1 

PARAMETER  PERTURBATION  METHODS 


(M 


The  adaptation  processes  in  these  systems  are  driven  by  the  error,  e, 
between  the  model  and  system.  This  error  is  used  to  form  a  performance 
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criterion,  which  for  this  example  was  integral  of  error  squared,  j  e  d  t, 
This  criteria  contains  a  component  of  the  frequency  at  which  a  parameter 
of  the  system  in  Fig*  1 (a)  is  being  perturbed.  The  amplitude  and 
phase  of  this  component  give  the  magnitude  and  sign  of  the  signal 
obtained  through  perturbation  and  integration.  This  signal  is  then 
fed  back  negatively  to  reduce  the  error.  The  adaptive  loop  mechan¬ 
ization  is  illustrated  in  Fig.  2. 


Fig.  2 

PARAMETER  PERTURBATION  ADAPTIVE  LOOP 


For  analog  simulation  electronic  multipliers  have  been  utilized  for 
parameter  correction  insertion.  In  operating  systems  servo-multi- 
pliers  could  be  utilized. 

The  system  of  Fig,  1(a)  suffers  from  the  disadvantage  that  the 
perturbation  appears  in  the  output.  The  system  of  Fig.  1(b)  avoids 
this  problem  through  model  parameter  perturbation.  This  second  method, 
however,  requires  that  the  model  be  nearly  the  same  form  as  the  system. 
Both  methods,  moreover,  require  phase  correcting  networks  to  correct 
the  phase  shift  between  parameter  perturbation  signal  in  the  model  and 
the  oscillator  signal.  Through  utilization  of  further  model  forms 
with  parameter  perturbation  the  authors  demonstrate  that  it  is  possible 
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to  prepare  a  configuration  which  is  both  system  and  signal  adaptive, 
using  multiple  models  with  perturbation  of  either  model  or  system 
parameters,  Rajaraman  in  a  later  article,  Ref,  4,  has  shown  the 
detailed  operation  of  the  signal  and  system  adaptive  process.  In 
summary,  this  signal  and  system  adaptive  process  requires  two  models 
and  two  adaptive  loops.  The  two  models  are  of  different  nature.  The 
first  is  an  ideal  version  of  the  process,  and  the  second  is  of  the 
same  order  as  the  process.  The  second  model  is  optimized  for  the 
input  by  the  first  adaptive  loop  operating  upon  the  difference 
between  the  output  of  the  two  models.  The  second  adaptive  loop 
causes  system  parameters  to  follow  the  optimized  second  model  par¬ 
ameters,  and  at  the  same  time  corrects  for  disturbances. 

Model  Output  Modulation 

Dorf  and  Byers  have  shown  the  feasibility  of  the  model  reference 
self  adaptive  system  using  output  modulation  techniques.  The  analysis 
of  the  techniques  and  a  discussion  of  experimental  results  are  found 
in  Ref,  5.  They  are  briefly  summarized  here.  Fig.  3  presents  the 
block  diagram  of  the  system. 

Fig ,  3 

MODEL  MODULATION  ADAPTIVE  LOOP 


=  10- 


In  this  adaptive  loop,  the  output  error  is: 


E 


(t) 


Cs  (t)  Cm  (t) 


11  +  A  sin  cd 
m 


2.1 


For  this  example  error  squared  has  been  used  as  the  criterion*  Thus, 
after  being  shaped  and  demodulated  a  phase-sensitive  d.c*  signal  is 
obtained.  When  this  signal  is  subjected  to  integration  and  filtering, 
the  open  loop  output  of  the  adaptive  circuit  becomes: 


2 

A  C 


m 


(C 


m 


V 


2.2 


The  adaptive  circuit  loop  was  closed  by  this  signal  being  inserted 
into  the  system  as  a  driving  function,  to  correct  a  system  parameter* 
The  result  of  this  multiloop  circuit  is  a  complex,  non-linear  system 
which  has  adaptive  characteristics  to  maintain  an  approximately  in- 
variant  system  transfer  function*  In  functional  notation  the  system 
output  might  be  described  as  a  function  of  the  following  variables: 

c(t)  “  f  (KAd>  V  |R|’  V  T>  M) 

K^d  ~  adaptive  loop  gain 

=  frequency  of  input  signal 

R  «  magnitude  of  input  signal 

=  filter  network  time  constants 
T  -  system  transfer  function,  uncorrected 

M  55  model  transfer  function 


The  adaptive  loop  filter  was  designed  to  filter  out  the  normal  system 
operating  signals*  The  speed  of  adaptation  depends  partially  on  the 
filter  transfer  function.  The  gain  and  lag  of  the  filter  influence 
the  magnitude  of  the  adaptive  response  and  the  input  signal  contami- 
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nation  One  important  feature  of  this  method  is  that  an  absolutely 
constant  frequency  oscillator  is  not  required  The  operation  of  the 


adaptive  loop  is  not  dependent  upon  frequency  stability  of  the  modu¬ 
lating  frequency.  Furthermore,  as  there  are  no  significant  time 
constants  between  modulation  and  demodulation,  the  frequency  of 
modulation  is  limited  only  in  bounds  at  the  lower  frequency  by  the 
necessity  of  being  an  order  above  the  input  signal  frequency.  Dorf 
and  Byers  further  investigated  a  two  parameter  correction  situation, 
driving  one  adaptive  loop  with  output  and  the  second  adaptive  loop 
with  the  derivative  of  output.  In  their  example  the  system  open 
loop  transfer  function  was: 


2  3 


The  model  open  loop  transfer  function  was  similar.  The  first  adaptive 


loop  was  driven  by  and  to  correct  (K^  +  A  K)  ,  the  open  loop  gain. 

The  second  adaptive  loop  was  driven  by  C  and  C  to  control  the  inner 

m  s 


loop  damping  (2  £  +  A  2  £  o^) .  The  authors  concluded  that  a  two 

parameter  self  adaptive  system  of  this  kind  was  feasible.  A  major 
restriction  was  found  to  exist  if  the  two  parameters  being  adapted 
were  not  independent,  and  the  adaptation  would  be  successful  only  if 
the  range  of  required  adaptation  was  restricted. 

Application  of  Model  Modulation 

To  verify  the  results  of  model  modulation,  and  to  obtain  exper¬ 
ience  with  the  characteristics  of  the  method,  the  method  was  simulated 
on  a  Donner  3100  analog  computer.  Initial  investigation  concerned 
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the  use  of  a  second  order  system  with  one  variable  parameter  being 
adapted  to  follow  a  second  order  model.  The  model  and  the  system 
were  ordered  for  an  open  loop  transfer  function: 


Gcs^  - 


-ort 


loo 


S  (s  tZ  S  (  S  4-  i  0s) 

The  parameter  drift  then  is  represented  as 


2.4 


G(V>  - 


“Ur 


iV 


s(s  +  +  »  ) 

where  n  was  a  step  parameter  drift.  The  input  was  also  of  various 
wave  forms  of  frequency  less  than  the  bandwidth.  The  adaptive  system 
open  loop  is  represented  as: 


2.5 


G  -- 


s(s  +  +  71  -  -™c) 

where  m^  is  the  correction  derived  from  the  output  driven  adaptive 
loop.  This  system  might  be  considered  to  simulate  a  second  order 
position  servo  as  displayed  in  Fig.  4. 

Fig.  4 

SECOND  ORDER  SERVOMECHANISM 


2.6 


Constant  Armature  Current  D.  C.  Motor  Drive 
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Thus  the  variable  drift  simulated  a  system  with  varying  or  unknown 
triction(f),  a  system  with  constant  gain  but  varying  velocity  error 
coefficient*  The  correction  would  be  inserted  by  varying  tachometer 
gain  feedback  magnitude . 

For  operational  simplicity  the  error  criterion  was  altered  from 
the  integral  of  error  squared,  ISE,  to  the  integral  of  absolute 
error  or  IAE.  Additionally,  the  modulation  and  demodulation  was 
accomplished  through  use  of  electronic  multipliers  *  Incidental  noise 
problems  with  these  multipliers  was  such  that  the  pure  integration 
followed  by  a  filter  was  removed  from  the  adaptive  loop,  and  instead 
a  double  filter  network  substituted.  The  alteration  is  shown  symbol¬ 
ically  in  Fig*  5. 


Fig.  5 

MODIFIED  MODEL  MODULATION  LOOP 


Satisfactory  results  were  achieved  with  this  circuit  with  the  filter 
break  points  a  =  . 1 ,  and  b  =  .1.  The  full  computer  diagram  is  shown 
in  Fig.  6  The  value  of  the  gain  of  the  adaptive  loop  for  the  circuit 
of  Fig  6  was  3.  Fig.  7  presents  the  output  data  for  this  same  circuit, 
illustrating  adaptation  for  various  parameter  drift  conditions 
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Additional  investigation  was  made  of  the  control  of  a  third  order 
system  with  one  variable  parameter  by  a  second  order  model,.  The  third 
order  system  investigated  open  loop  transfer  function  was 


G(s) 


K  100 

____£_ _ _ _ _ _ _ _ _ 

s  (s  +  K  )  (s  +  10  +  n  -  m  ) 

p  v  c 


2.7 


The  model  utilized  was  that  of  the  previous  example.  The  analog 

computer  setup  was  that  of  the  transfer  function  format  similar  to 

K 

Fig.  6,  with  an  additional  function  p  installed.  The  adaptive 

s+K 

P 

loop  was  unchanged.  This  system  then  simulated  a  third  order  control 
system  as  illustrated  in  Fig*  8. 
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Fig.  8 

THIRD  ORDER  SYSTEM 


For  this  situation  the  variable  parameter  simulated  would  be  p, 
with  the  correction  being  inserted  through  a  small  servo  varying  m^. 

The  adaptation  was  satisfactory  for  =  40 ,  when  p  =  10.  Adapt¬ 
ability  became  difficult  for  <  40  for  large  parameter  shift.  Fig.  9 
presents  the  output  data  for  this  circuit  for  a  representative  =  40 , 
under  various  parameter  drift  conditions. 

As  a  final  illustration  of  the  model  modulation  method,  cancel¬ 
lation  techniques  were  used  for  system  control.  The  process  is  best 
illustrated  in  Fig.  10. 
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Third  Order  Svsteoa  Adapting  to  Second  Order  Model 
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Fig.  10 

CONTROL  SYSTEM  CASCADE  COMPENSATION  ADAPTATION 


This  technique  is  most  useful  when  the  operating  system  pole  is  not 
amenable  to  correction  because  of  the  physical  characteristics  of  the 
system.  In  the  operation  illustrated  by  Fig.  10  zero  drift  caused 
compensator  cancellation.  When  drift  was  introduced  the  correction 
forced  the  zero  of  the  compensator  to  cancel  the  system  pole,  and  the 
desired  system  operation  was  restored.  The  analog  computer  operating 
schematic  is  presented  in  Fige  11.  For  operational  simplicity  the 
drift  was  inserted  as  a  step  change  in  the  system  pole.  The  computer 
output  data  is  presented  in  Fig.  12,  where  the  successful  adaptation 
of  the  system  is  well  illustrated. 
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Fig.  11 

ANALOG  COMPUTER  SIMULATION  OF  COMPENSATION  ADAPTATION 


The  results  of  these  simulations  of  the  adaptive  system  were 
considered  to  satisfactorily  demonstrate  the  feasibility  of  the 
method  and  to  serve  for  participant  training.  The  analytical  proof 
of  the  correctness  of  the  alterations  to  the  initial  formulation, 
IAE  for  ISE  as  criterion,  and  filter  instead  of  pure  integration, 
have  been  neglected.  The  demonstrated  accuracy  of  the  adaptation 
was  felt  to  be  sufficient  proof  of  the  method. 
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Chapter  III 


IDENTIFICATION  THROUGH  MODEL  MODULATION 

The  methods  of  successful  adaptive  control  through  model  modu¬ 
lation  were  directly  applicable  to  the  identification  of  operating 
systems.  The  requirements  initially  set  were  that  the  identification 
should  proceed  utilizing  only  the  operating  system  input  and  output 
without  disturbing  the  system  operation  lo  utilize  the  model  rnodu- 
lation  method,  the  model  would  be  placed  in  parallel  operation  with 
the  system.  The  adaptive  loop  output  would  be  applied  to  the  model 
parameters  to  cause  the  model  to  adapt  to  the  system  Thus,  knowledge 
of  the  unadapted  model  governing  equations  combined  with  the  magnitude 
of  the  adaptive  loop  output  would  serve  to  define  the  system. 

This  method  is  shown  in  block  diagram  form  in  Fig.  13 

Fig.  13 

MODEL  MODULATION  WITH  MODEL  ADAPTATION 


This  method  of  system  identification  does  not  relieve  the  engineer 
of  the  necessity  of  system  analysis.  The  form  of  the  transfer  function, 
the  order  of  the  operating  equation,  and  the  general  magnitude  of 
the  poles  and  zeros  of  the  system  must  be  known  in  order  to  synthesize 
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the  model  Ideally,  one  or  two  critical  parameters  whose  magnitude 
may  not  be  susceptible  to  analysis  may  be  estimated,  then  measured 
or  identified  by  the  model  modulation.  These  limitations  should  be 
expected  from  any  adaptive  system.  The  knowledge  required  to  build 
an  adaptive  loop  which  would  assume  no  knowledge  of  the  system,  and 
then  synthesize  perfect  identification,  with  sufficient  economy  to 
be  practical,  is  beyond  the  state  of  the  art 

The  standard  second  order  servo  or  feedback  system  equation  is 
formulated  in  transfer  function  terms  as; 


Cm  _  _ 

t 


For  this  system  the  standard  open  loop  transfer  function  is: 


3. 1 


z  £  (S3  =  _ _ _  3. 

L(s\  5  (s  *  2 

The  response  of  the  second  order  system  is  completely  described  by 

the  two  variables,  zeta  and  uo  A  higher  order  system,  or  any 

n 

system  with  closed  loop  transfer  function  zeros,  naturally,  has 
additional  significant  variables.  If,  however,  it  is  considered  that 
there  are  two  variables  of  primary  interest,  the  higher  order  open  loop 
transfer  function  may  be  written  as: 

_  ns  _  -^2  k£  3. 

4  s(stZ%)  ft  (nf <  ) 

v  K~i  1  K 


where  u>  and  zeta  are  these  variables, 
n 
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In  the  development  of  model,  modulation  system  identification  method 
two  major  assumptions  were  made.  The  first  of  these  is, 

1 

That  in  the  operating  system  to  be  identified  there  will  be  but 
two  unknown  parameters.  Thus  the  system  must  be  second  order, 
a  completely  dominant  second  order  approximation  oC  a  higher 
order  system,  or  a  higher  order  system  which  has  the  ether 
parameters  fully  identified  If  the  latter  is  the  case,  the 
format  of  the  preceding  equation  will  be  used  to  describe  the 
system. 

This  first  assumption  limited  the  parameter  adaptation  problem  to  the 
two  variables  zeta  and  ^ 

Examination  of  any  physical  system  reveals  that  zeta  and  may 

or  may  not  be  functionally  independent  The  simple  second  order  system- 
detailed  in  Fig.  5  serves  as  an  adequate  example.  The  closed  loop 
transfer  function  is: 

CH)  -  rj 

.  s6'  t  s  f  5-  +  t'l  +  i<Jku 

1  J  rj  >  rj 

In  this  system  variation  of  the  friction(f)  will  yield  an  independent 

zeta  and  id  „  Variation  of  all  other  factors  will  result  in  zeta  and 
n 

U)  being  different  functions  of  the  same  variable,  and  thus  not  in- 
n 

dependent  Thus  zeta  and  u>  may  be  independent  or  both  different 

n 

functions  of  the  same  variable. 
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The  adaptive  method  that  has  been  discussed  to  this  point  has 

dealt  with  but  one  variable-  Since  the  interdependence  of  zeta  and  U) 

n 

is  at  the  most  one  of  their  being  two  different  functions  of  one  factor, 
it  is  reasonable  to  expect  that  two  adaptive  loops  will  be  required 

The  second  major  assumption  required  depends  heuri s ti ca 1 ly  on  this 

discussion  of  the  functional  relationship  of  zeta  and  od  c  This 

n 

assumption  is : 

II. 

That  model  may  be  adapted  to  the  system  through  two  separate 
adaptive  loops.  The  first  of  the  loops  will  operate  upon  the 
system  and  model  output  to  develop  the  first  corrective  factor 
The  second  adaptive  loop  will  operate  upon  the  derivative  with 
respect  to  time  of  the  system  and  model  output  to  develop  the 
second  corrective  factor. 

The  remainder  of  this  report  is  devoted  basically  to  showing  that  these 
two  assumptions  are  justified 


Proceeding  on  these  assumptions  the  investigation  utilized  models 
of  the  following  general  form: 


k6*&-  «c)  (St 

s(s  ttSrt,,  +-"d)  Tr  ) 


3.5 


A  schematic  of  the  operation  of  these  loops  as  tentatively  applied  to 
system  identification  is  illustrated  in  Fig.  14. 
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Fig.  14 

SCHEMATIC  OF  TWO  PARAMETER  CORRECTION 


With  this  model  the  characteristic  equation  of  the  operating 
system  could  then  be  easily  detailed.  For  the  model  with  the  form 
of  equation  3.5  and  unity  feedback,  the  closed  loop  transfer  function 
becomes : 


This  function  may  then  be  used  to  accurately  detail  the  system  differ* 
ential  equations. 


Only  limited  single  parameter,  single  adaptive  loop  correction 
investigation  was  made.  Primary  investigation  concerned  the  two  par 
ameters  of  the  previous  equation  as  in  Fig,  14 

Fig.  14  also  illustrates  the  visualized  mechanization  of  the  model 
modulation  identification  system  Potentially  the  model  and  the  adaptive 
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loops  would  be  constructed  through  use  of  standard  operational  amplifier 
networks.  Thus  any  moderate  sized  analog  computer  operating  in  con¬ 
junction  with  the  scaled  input  and  output  of  a  physical  system  would 
serve  as  an  identifier. 
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Chapter  IV 


ANALOG  COMPUTER  INVESTIGATION 


Model  Modulation  1  den t i f i ca t ion  o t  an  E qua t ion _o|_ Known  Order 

The  primary  investigation  was  centered  around  use  of  the  analog 
computer,  The  investigation  concerned  use  of  two  undetermined  par 
ameters  in  the  system.  The  simulated  control  system  transfer  function 
was  of  the  form  shown  in  formula  3,3  The  model  transfer  function  as 
of  the  form  shown  in  formula  34,  In  both  cases  the  actual  computer 
arrangement  was  that  of  the  transfer  function  format  The  system 
values  of  u)  and  zeta  were  set  up  through  switching  circuits  so  that 
three  different  values  of  each  could  be  used  directly  The  adaptive 
loops  used  on  the  computer  were  those  shown  in  Fig  6  In  general, 
both  filters  in  the  adaptive  loops  were  set  for  the  break  at 
omega  ®  A  radian.  Fig.  15  presents  the  full  analog  computer  diagram 
for  a  simulated  third  order  system  with  a  model  tracking  two  unknown 
parameters o  This  diagram  conforms  to  the  format  of  formula  3,3  and  3  5. 
The  initial  investigation  with  this  computer  a rrangement  concerned  a 
control  system  with  simulated  open  loop  transfer  function  of: 


G(s)  - 


10  (100)  K 
s  (s  +  p)  (s  +  100) 


5  £  K  £  15 
5  *  p  *  15 


4, 1 


and  a  model  transfer  function: 
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Flrure  15 

Analog  Computer  Diagram  for  Two  Parameter 
Model  Adepts t ion 


The  input  was  primarily  deterministic.  It  consisted  of  a  square 


wave  with  a  frequency  which  was  less  than  the  bandwidth  Tests  were 

also  carried  out  however,  using  a  random  input.  This  random  input 

was  the  output  of  a  Gaussian  "white"  noise  generator  with  an  18  db 

cut  off  at  approximately  the  system  bandwidth.  At  the  completion  of 

each  run,  the  recorded  m^  and  m.  were  used  in  conjunction  with  the 

known  model  characteristics  to  derive  the  closed  loop  characteristic 

equation  of  the  simulated  system.  This  allowed  comparison  of  actual 

system  zeta  and  Cl)  with  those  derived  from  that  identified  as  the 
n 

system.  A  non-dimensional  corrective  parameter  was  examined  and 
used  as  a  partial  criterion  of  identification  This  corrective 
parameter  a^  is  best  understood  through  examination  of  Fig.  16. 


Fig .  16 

CORRECTIVE  PARAMETER 
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In  this  figure,  (T)  is  the  dominant  root  location  of  the  simulated 
control  system,  (2)  is  the  root  location  of  the  uncorrected  model, 
and  (3)  is  the  position  of  the  corrected  model  root.  The  arithmetic 
mean  if  and  the  standard  deviation,  o  of  a^,  were  computed  for  each 
series  of  tests. 

-  36- 


The  data  for  the  first  series  of  tests  is  shown  in  Table  I  In 


this  table  the  variance  of  the  tracked  parameters  is  shown.  The 

roots,  zeta  and  U)  have  been  shown  for  only  the  two  extreme  cases 
n 

in  each  test,  that  of  the  maximum  and  minimum  zeta.  The  correction 
parameter  a^  and  the  sigma  have  been  calculated  from  the  assembled  data 
of  the  eight  situations  for  each  test  where  the  parameters  do  not 
match.  Fig.  17  shows  a  sample  of  the  analog  computer  output.  This 
data  is  that  of  the  first  and  last  cases  tabulated  in  Table  1 
Figure  17D  will  be  used  to  illustrate  the  method  of  computing  the 
characteristic  equation  roots  of  the  operating  control  system  For 
this  example,  the  model  closed  loop  transfer  function  was 

10.7  (s+14)  (10-m  ) 

R  (s)m  =  s  (s+15)  (s+l(Hnu)  +  10*7  (s+14)  (10-m  ) 

From  figure  17D  the  values  of  the  correction  factors  may  be  seen  to 

be  m  =  1.5  and  m*  -  1.080.  The  characteristic  equation  is  then 
c  c 

s3  +  26.08  S2  +  289. 1  s  + 1724.  44 

The  roots  of  this  equation  were  obtained  through  digital  computer 
service.  These  roots  are: 

-5.894  +  j  9.268 

-14.29 

Fig.  18  shows  a  plot  of  the  characteristic  equation  dominant  root 
positions  for  the  actual  system  and  the  identified  system  detailed 
in  Fig.  17. 
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Equation  Order  Approximated 


During  the  process  of  analyzing  an  operating  system,  it  is 
inevitable  that  approximations  would  be  made.  These  approximations 
will  generally  mask  or  ignore  the  very  small  time  constants,  which 
in  the  s  plane  are  the  poles  of  large  magnitude.  In  trying  to  use 
the  model  modulation  method  to  identify  an  operating  system  it  is 
then  inevitable  that  the  model  and  the  system  would  be  of  different 
order.  For  a  step  input,  the  effect  of  these  poles  would  be  greatest 
during  the  initial  portion  of  the  transient  response,  and  the  effect 
would  become  inconsequential  in  the  later  portions  of  the  response. 

It  was  hypothesized  that  if  the  uncertain  parameters  of  the  operating 
system  concerned  dominant  roots,  satisfactory  identification  would 
take  place.  The  restriction  to  this  hypothesis  was  two  fold.  It 
should  be  expected  that  if  the  input  was  of  a  random  noise  nature 
with  no  stead  state,  no  identification  would  be  possible  unless,  as 
in  the  first  test  of  Table  I,  the  equation  order  and  the  additional 
parameters  are  known.  Additionally,  with  deterministic  inputs, 
identification  could  take  place  only  when  these  unknown  poles  were 
of  sufficient  magnitude  to  be  unimportant  with  respect  to  the  dominant 
response . 

To  examine  this  restriction  the  tests  detailed  in  Table  II  were 
hold.  Since  the  systems  to  be  tested  were  set  up  in  transfer  function 
format  on  the  analog  computer,  each  pole  was  a  new  block  set  Into  the 
basic  arrangement  of  Fig.  14.  The  information  in  Table  II  is  of  the 
same  nature  as  in  Table  I,  but  as  all  nine  system  configurations  differed 
from  those  of  the  model,  all  nine  situations  were  used  to  calculate 
the  a  and  o. 
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The  results  of  these  tests  were  as  expected  The  random  input 
revealed,  as  expected,  lack  of  identification.  Under  deterministic 
inputs,  when  the  postulated  restrictions  were  held,  satisfactory 
identification  took  place.  One  interesting  point  revealed  in  this 
investigation  was  that  an  adequate  test  for  satisfactory  identification 
would  be  made  if  the  frequency  of  the  input  could  be  varied.  Naturally 
this  violates  the  basic  restriction  that  led  to  the  development  of 
this  method,  i.e.,  no  violation  of  the  natural  operation  of  the  system, 
but  is  a  point  of  interest.  For  this  process  the  frequency  was  varied 
between  a  lower  frequency  with  a  period  equal  to  the  second  order 
time  constant  and  an  upper  frequency  equal  to  the  bandwidth.  When 
the  adaptive  loop  outputs  remained  relatively  constant  through  this 
frequency  range  it  was  found  that  these  outputs  served  to  adequately 
identify  the  dominant  roots  when  the  characteristic  roots  were  found. 
If,  however,  the  identification  varied  through  this  range,  the  effect 
of  the  ignored  poles  precluded  satisfactory  identification. 

Figure  19  presents  illustrative  data  from  the  first  tablulated 
case  of  Table  II  for  this  series  of  tests.  Fig.  20  shows  the  dominant 
roots  of  the  system  for  each  of  the  tests  illustrated  in  Figc  19, 
with  the  dominant  roots  as  identified  from  the  data. 

Examination  of  the  data  reveals  that  for  both  of  these  types  of 
equation  situations  good  dominant  root  identification  has  been  achieved 
It  was  found  that  these  dominant  roots  could  be  identified  to  within 
an  average  error  of  2 . 27>  of  the  undamped  natural  frequency  and  to 
within  an  average  error  of  4.767>  of  the  complex  phase  angle.  When 
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this  method  is  used  to  identify  the  characteristics  of  an  actual  linear 
control  system,  it  is  highly  probable  that  greater  accuracy  may  be 
attained  since  the  model  parameters  may  be  varied  to  obtain  error 
bracketing.  This  process  is  illustrated  in  Chapter  VI  * 
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Chapter  V 


DIGITAL  COMPUTER  SIMULATION 

Direct  utilization  of  a  digital  computer  for  identification  of 
an  operating  control  system  was  not  considered  1 t  is  possible, 
however,  that  operating  systems  might  be  monitored  or  controlled  by 
digital  techniques.  In  this  case  the  output  would  be  presented  in 
the  form  of  discrete  numerical  data.  Ideally  this  data  could  be  used 
as  input  to  a  digita 1- to-analog  converter,  and  the  analog  techniques 
previously  described  could  be  used  for  operating  system  identification. 

In  the  absence  of  this  equipment,  digital  simulation  of  the  control 
system,  model,  and  adaptive  loops  was  investigated  to  illustrate  the 
techniques  and  the  potential  applicability  of  model  modulation  techniques 
through  numerical  techniques.  This  investigation  was  conducted  through 
the  use  of  the  CDC  1604  Digital  Computer.  The  program  code  was  written 
in  Fortran.  A  standard  Runge-Kutta  numerical  integration  process  was 
used  to  solve  the  differential  equation  which  governed  the  operation 
of  the  full  model,  system  and  adaptive  loop.  The  differential  equations 
were  formulated  in  a  transfer  function  format 

Three  basic  configurations  of  model  adaptation  were  tested  through 
this  digital  computer  simulation  For  these  configurations  the  model 
was  used  to  identify  the  operating  control  system  characteristic 
equation  roots  under  the  following  conditions 

1.  Second  order  control  system  with  undetermined  pole. 

2.  Second  order  model,  third  order  control  system  with  uncertain 
gain  and  uncertain  dominant  pole. 

3.  Third  order  system  and  model  with  two  uncertain  poles. 
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The  actual  investigation  was  limited  to  achieving  a  successful  solution 
for  one  parameter  set  for  each  of  these  configurations.  The  excessive 
time  used  by  integration  techniques  on  the  digital  computer  was  of 
course  a  limiting  factor.  The  first  of  these  configurations  will  be 
by-passed  in  this  discussion  since  it  can  be  considered  to  be  a  sub¬ 
system  of  the  second  configuration. 


The  model  for  the  second  configuration  had  the  following  open 
loop  transfer  function: 


10  (10  -  m  ) 

c 

G(s')  =  s  (s  +  10  +  m.) 

c 

and  a  closed  loop  transfer  function  of 

r  10  (10  -  m  ) 

R  (S;m  "  s  (s  +  10  +  m.)  +  10  "(10  -  m  ) 

c  c 


5.  1 


5.2 


The  control  system  open  loop  transfer  functions  were 


G  (s'! 


/£  ooo 

S  (s  +  loo)(s  +  8) 


5.  3 


and 


(s')  ^  _  5  4 

S  (  S-J-  ioo)(  S  +  1 2  ^ 

Figure  21  presents  the  transfer  function  format  which  was  used  to 

formulate  the  differential  equation  for  the  FORTRAN  program.  This 

format  allowed  formulation  of  the  operation  by  a  set  of  quasi-state 

variables  of  the  form: 


x 

n 


f  (x  >  x  ,  x  ,  000  x  ,  R) 
n  1  2  J  n 


5  5 
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The  solution  of  the  equations  was  programmed  to  be  presented  in  both 


numerical  and  graphical  format.  Fig.  22  presents  a  sample  of  the 
input  program  Fig.  23  presents  a  sample  of  the  numerical  data  output, 
and  Fig.  24  presents  the  graphical  data  output  for  the  control  system 
open  loop  transfer  function  given  by  equation  5  3. 

The  graph  of  Fig.  24  presents  the  corrective  parameters  as 
"corrected  pole"  =  (10  +  m.)  and  "corrected  gain"  =  (10  -  m^) 

Table  III  shows  the  identification  of  the  characteristic  equation 
roots  for  these  two  situations. 


Table  IIIA 


DIGITAL  COMPUTER  CHARACTERISTIC  EQUATION  ROOTS 


Control  System 


Model  as  Corrected 


Equation  5 . 3 


3.36  +  j  10. 35 
101.27 


3.07  +  j  997 
101.15 


Equation  5.4 


5.55  +  j  6.96 
100.9 


4.88  +  j  7.77 
100.93 
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Figure  2.2.  a. 


Digital  loaputer  Program 


PROGRAM  ADAPT  12 

C  N=NO .OF  EONS.,  NP=NO.  OF  INCREMENTS  BETWEEN  PRINT  OUTS,  TO=INITIAL 
C  TIME,  TF  =  FINAL  TIME,  CT  =  T IM  E  INCREMENT,  TN=  T I  ME  FIRST  PRINTED  DATA, 

C  X ( J ) *  I N  I  T I AL  VALUES  OF  THE  VARIABLES  X ( J ) 

C  NG=NO.  OF  INCREMENTS  BETWEEN  GRAPH  POINTS 

OD I  MENS  I  ON  X ( 30) , Xl ( 900)  ,X2 (900 ) ,X3( 900) ,X4( 900) , X5( 900 ) , X6 ( 900 ) , T I 
1  ME (900),  H ( 20 ) 

COM  MCN  R,  E,  T  2 ,  H 
100  FORMAT  (31  1C/(PF10.0) ) 

102  FORMAT  (  12  /  (  I2.F2C.0  )  ) 

200  FORMAT  (F8.3,1CF11.5) 

2010F0RMAT  ( 99H  TIME  SYS. PCS.  MOD.POS.  POS.DIFF.  SYS.VEL, 


GAIN  COR.  POLE  COR. 


//) 


2H).  /  ) 

,  27H,  POLE  CORR. 


1  MOD  .  VEL .  VEL.DIFF. 

202  FORMAT  (/) 

203  FORMAT  ( 1 H 1 > 

204  FORMAT  (//) 

205  FORMAT  ('9H  SYSTEM^  F9.  1,  1  3H/ S ( S  + 1 00 ) ( S  +  F6.2, 

2060F0RMAT  (,32H  GAIN  CORR  FILTER  CORNER  AT  W  =  £6.2 

1  FILTER  .£T  W  =  F6.2  /  )  | 

207  FORMAT  f25H  MODEL  =  1  CO . / S ( S+ 1 0 ) .  v/  ) 

2080FORMAT  (37H  ADAPTIVE  SYSTEM  GAINS,  OUTPUT  (S)  =  F6.2,  20H, FIRST  DE 
1 R I  V .  ( D)  =  F6.2  /) 

210  FORMAT  ( 20H  COEFFICIENT  NUMBER  12,  2H  =  F12.4,  1H.  ) 

READ  100,  N,  NP,  NG ,  TO,  TF,  DT ,  TN ,  (X(J),  J=1,N) 

READ  102,  NF,  (  K,  (H(K)),  I  =  1 , NH  ) 

PRINT  204 

PRINT  208,  H  16),  H { 1 7) 

ASYSG  =  H ( 5  )  *  H  (  2  )  «  H{3 ) 

PRINT  205  ,  ASYSG,  H(4) 

,  H (  11),  H(  15  ) 


I  =  1  ,  NH 


1 


PRINT  206 
PRINT  207 
DO  6 

PRINT  202 
PRINT  210, 
PRINT  203 
PRINT  201 
T  =  TO 

T2  =  0.0 
NUMPTS  =  0 
MP T  S=  0 
DPOS  =  XU) 
DVEL  =  X ( 2 ) 


I  ,  h  (  I  ) 


X  (  4  ) 

X  {  5  ) 

AGACOR  =  H ( 7 )  -  X(7)*H( 16) 

APCLCOR  =  H ( 6  )  * (  1.  +  X  { 9  )  *H (  17)) 

IF  (T-TN)  36,  2,  2 
IF(MPTS)  3C »  35 , 30 
IF ( XMODF (MPTS, 50*NP ) )  31,33,31 
IF ( XMODF ( MPTS, 10*NP) )  32,34,32 
IF ( XMODF ( MPTS,  NP ) )  36,35,36 

33  PRINT  203 
PRINT  201 

34  PRINT  202 

35  PRINT  200,  T , X (  1  )  ,  X ( 4  ) , DPOS , X ( 2 ) , X ( 5 ) ,  DVEL,  AGACOR , APOLCOR 

36  MP  T  S=MP  TS  +  1 

37  I F ( XMODF (MPTS,  NG ) )  38,  39,  38 

39  NUMPTS  =  NUMPTS+1 

IF  (900  -  NUMPTS  )  38,  38,  40 

40  XI ( NUMPTS)  =  X ( 1 ) 


2 

30 

31 

32 


X2 ( NUM  PT  S ) 

X3(  NUMIPTS  ) 

X4 (NUMPTS ) 

X5( NUMPTS ) 

X6( NUMPTS) 
TIME  (NUMPTS) 


X  (  4  ) 
AGACOR 
APOLCOR 
DPOS 
DVEL 
=  T 


38  IF ( TF-T-DT )  10,20,20 

10  PRINT  204 


CALL 

GRAPH 

( 

NUMPTS 

» 

T  IME 

9 

X  3 

9 

8 

) 

CALL 

GRAPH 

( 

NUMPTS 

f 

T  IME 

9 

X  4 

9 

8 

) 

CALL 

GRAPH 

( 

NUMPTS 

9 

T  IME 

9 

X5 

9 

8 

) 

CALL 

graph 

( 

NUMPTS 

9 

T  IME 

9 

X  6 

9 

8 

) 

CALL 

GRAPH 

( 

NUMPTS 

9 

T  IME 

9 

XI 

9 

8 

) 

CALL 

GRAPH 

( 

NUMPTS 

9 

T  IME 

9 

X2 

9 

8 

) 

STOP 

20  CALL  RKUTTA (N »  T,  X,  DT  ) 
T2=  T  2  +  DT 


Fi.m  rfi  £.2.  b 


Digital  Computer  Program 
( con ' t ) 


C 

C 


r  C  (  14  >  ,  H  (  30  ) 


tNU 

SUBROUTINE  RKUTTA  (NtT,X,DT) 

DIMENSION  X ( 30  )  tAK(4,30)  ,  XDCTI30),  XC(30),I 
COMMON  R •  E  t  T2f  H 
C ( 1  )  =0.0 
C  I  2  )  =0.5 
C ( 3  )  =0 . 5 
C ( 4 )  =  1  .0 
DO  4  1=1,4 

TC=T  +  ri(  I  )  «  DT 
DO  2  J=T.N 

XC( J)=  XI  J)  +CII)*  AK I  I  - 1 , J  ) 

CALL  DERI VI TC ,XC,XDOT  ) 

CO  4  Jf 1 v N 

AK  (  ItJ)=  DT*  XCOT(J) 

DO  3  J  =  1  »  N 

X( J)=X( J)  +  ( AKC1 ,  J)  +  2.*AK(2, J)+2.»  AK I  3 , J ) + AK ( 4 , J ) ) /  6. 

RETURN 

END 

SUBROUTINE  CERIVIT.  X,  XDOT) 

DIMENSION  XDOTI30),  XI3C),  H  (  30 ) 

COMMON  R,  E,  T2,  H 
I  F ( 1 .0-  12)3,3,1 
R  =  H(  1  ) 

GO  TO  2 
R  =  -H( 1) 

I F { 2 . 0  -T2 )  4,4,2 
T2=  0.0 
GO  TO  1 

E  =  SINF(314.*T ) 

H ( 4  )  =  VARIABLE  POLE 
H ( 5  >  =  VARIABLE  GAIN 
XDOT (  1 )  =  H ( 5 ) *  X ( 2 ) 

XDOT ( 2 )  =  X  (  3  )  «H ( 3 )  -  X(2)*H(4) 

XDOT ( 3 )  =  H ( 2  )  * ( R  -  X ( 1  ) —  X ( 3 ) ) 

-  ‘  X ( 5 ) * ( H ( 7 )  -  X ( 7 ) *H ( 16)) 

H  (  6  )  * ( R  -  X ( 5 )  -  X  I  5 ) *X I  9 ) *H I  17)  -  X ( 4 ) ) 

H  (8  )  *E*A8SF(X ( 1 )-  X ( 4 )  -  E*X(4))  -  X(6)*H(9) 

HI  10  )  *X( 6)  -  X  t  7 ) *H ( 1 1 ) 

HI  12)*E*ABSF(  X ( 2 )  -  XI5)  -  X(5)*E)  -  X(8)*H(13) 
H  1 1 4  )  *X I  8  )  -  HI 15)*X(9) 


XDOT (4 ) 
XDOT I  5 ) 
XDOT I  6  ) 
XDOT ( 7 ) 
XDOT I  8 ) 
XDOT I  9 ) 
RETURN 
END 
END 


•-64- 


Figure  2-3  a. 

Digital  Computer  Numerical  Output 


ADAPTIVE  SYSTEM  GAINS,  OUTPUT  (S)  =  10. 00, FIRST  OERIV.  (D)  = 

SYSTEM*  12000. 0/S(S+lC0) { S+  8.00). 

GAIN  CORR  FILTER  CORNER  AT  W  =  .10,  POLE  CORR.  FILTER  AT  W 


MODEL  =  100. 

/SI S+  10  )  . 

COEFFICIENT 

f^flBER 

NUMBER 

1 

= 

10.0000. 

COEFFICIENT 

2 

= 

100.0000. 

COEFFICIENT 

NUM8ER 

3 

= 

1 0.0000. 

COEFFICIENT 

NUMBER 

4 

= 

8 .0000. 

COEFFICIENT 

NUMBER 

5 

S 

1 2.0000. 

COEFFICIENT 

NUM8ER 

6 

= 

1 o.cooo. 

COEFFICIENT 

NUMBER 

7 

= 

10.0000. 

COEFFICIENT 

NUMBER 

8 

.5000. 

COEFFICIENT 

NUMBER 

S 

.5000. 

COEFFICIENT 

NUMBER 

1C 

= 

. 1000. 

COEFFICIENT 

NUMBER 

1  1 

= 

.1000. 

COEFFICIENT 

NUM8ER 

12 

= 

.5000. 

COEFFICIENT 

NUMBER 

13 

S 

. 5000. 

COEFFICIENT 

NUMBER 

14 

= 

.  1000. 

COEFFICIENT 

NUMBER 

15 

= 

.  1000. 

COEFFICIENT 

NUMBER 

16 

S 

1 o.cooo. 

COEFFICIENT 

NUMBER 

17 

= 

1  0 «  0  0  0  0  © 

10.00 

.]0 
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Digital  Tomoutpr  Mu  Output 
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Digital  Omputer  Graphical  Output  fir  Two  Parameter  Identification 


The  third  configuration  examined  was  that  for  the  following  model 


open  loop  transfer  function: 


G(s) 


10,000 _ 

s  (s  +  100  +m  )  (s  +  10  +  m.) 
\  c'  x  c 


5  6 


with  control  system  open  loop  transfer  functions  of: 


G(s) 


10,000 _ 

s  (s  +  150)  (s  +  15) 


5.7 


and 


G(s) 


10,000 _ 

s  (s  +  75)  (s  +  7.5) 


5.8 


The  equations  were  formulated  in  the  manner  of  Fig  20,  except 
that  the  adaptive  loops  each  utilized  one  filter  and  one  integrator 
instead  of  the  two  filters  previously  detailed.  Satisfactory  results 
were  obtained.  Fig.  25  presents  sample  results  from  the  investigation 
of  equation  5.7,  and  Table  IV  presents  the  root  identification 
achieved. 


Table  IIIB 

DIGITAL  COMPUTATION  CHARACTERISTIC  EQUATION  ROOTS 
TWO  UNCERTAIN  POLES 


Equation  5.  7 


Control  System 

-  7.25  +  j  3. 72 

-  150.49 


Model  as  Corrected 

-7.87  +  j  2. 1 
-150.49 


Equation  5.8 


-  2.81  +  j  11.05 

-  76.87 


-3.45  +  j  11.05 
-  74.56 
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Digital  Couputer  Graphics!  Output  for  Integration  in  Ada-ptive  Loops 


Chapter  VI 


ACTUAL  SERVOMECHANISM  IDENTIFICATION 


Intro due t ion 

The  previous  discussion  and  the  data  in  Table  1  and  Table  II 
indicate  that  the  identified  location  of  the  roots  of  an  operating 
system  is  the  function  of  many  variables.  The  type  of  model  and 
the  proximity  of  the  model  equation  roots  to  the  roots  of  the  equation 
of  the  operating  system  are  obviously  factors  of  importance  in  this 
identification.  With  this  point  in  mind,  it  is  suggested  that  the 
following  steps  be  followed  in  the  process  of  identifying  an  actual 
system. 

1.  Analyze  system  operating  factors  to  enable  determination  of 
open  and  closed  loop  transfer  functions 

2.  Determine  which  values  in  the  open  loop  transfer  function 
are  questionable. 

3.  Utilize  operational  amplifiers  or  an  analog  computer  to 
prepare  the  dynamic  model  and  the  adaptive  loops 

4  Scale  the  model  input  and  outputs  to  the  estimated 
corresponding  values  of  the  system. 

5o  Set  indicated  adaptive  loop  polarities  to  adapt  the  model 
to  the  undetermined  or  questionable  values  in  the  system 

6.  Adjust  gain  of  adaptive  loop  to  obtain  satisfactory  operation 

7 :  Calculate  the  transfer  function  of  loop  corrected  model, 

scale  a  new  model  to  these  functions  and  repeat  operation 
to  generate  better  conformance  or  to  confirm  values  pre¬ 
viously  attained 
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A  position  servomechanism  was  used  to  illustrate  the  applicability 
of  the  model  modulation  techniques  to  an  actual  operating  system. 

The  servo  used  was  one  developed  for  use  in  instructional  laboratories 
by  Professor  J.  R.  Ward  of  the  U.  S .  Naval  Postgraduate  School,  and 
is  partially  described  in  Ref.  6.  The  full  description  and  details 
of  the  operation  of  this  system  are  on  file  at  this  school. 

Fig.  26  presents  a  block  diagram  of  the  system.  This  repre¬ 
sentation  shows  that  the  servomechanism  may  be  considered  as  a  gear- 
head  system.  Thus  the  use  of  this  identification  method  conforms  to 
the  stated  requirement  of  parallel  model  operation. 

Fig  26 

POSITION  SERVOMECHANISM 


Analysis  of  this  mechanism  resulted  in  the  transfer  function  diagram 
of  Fig  27. 
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Fig.  27 

POSITION  SERVOMECHANISM  BLOCK  DIAGRAM 


The  approximate  parameter  values  in  the  system  were  determined.  The 
magnitude  of  inductance  (L)  was  such  that  the  approximation  L  0 
was  made.  The  open  loop  transfer  function  was  then • 


£ 


G(s)  - 


RJ 


yn 


*( 


S  4- 


Rf  t  K>*  ) 
RJ  / 


)< 

J=>) 


6  1 


The  value  of  friction(f)  was  difficult  to  evaluate  because  of  the 
expected  non-linearity  of  friction  in  a  d.c*  motor.  A  nominal  value 
was  chosen  to  represent  the  viscous  friction,  while  coulomb  friction 
effects  were  ignored.  It  was  further  assumed  that  the  open  loop  gain 
was  questionable.  As  a  first  estimate  it  was  assumed  that  the  open 
loop  transfer  function  was: 


r  -  k  ^  400 

*  sTsTpT  s (s  +  8) 

Fig.  28  presents  the  control  sys tern-model -adaptive  loop  block 
diagram . 
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Fig.  28 

SERVOMECHANISM  IDENTIFICATION  BLOCK  DIAGRAM 


The  actual  model  and  adaptive  loop  analog  computer  schematic  was  that 
detailed  in  Fig.  15,  while  the  scaling  factors  were  achieved  through 
cascaded  operational  amplifiers  to  assure  correct  polarity.  The 
servomechanism  and  the  model  were  driven  by  a  square  wave  of  one  cycle 
per  second.  Table  IV  presents  the  results  of  four  tests  made  first 
to  bracket,  then  to  closely  approximate  the  open  loop  transfer 
function.  Fig.  29  presents  the  results  of  the  final  run  listed  in 
Table  IV.  One  the  basis  of  these  results  the  transfer  function  of 
this  servomechanism  was  estimated  to  be: 


e0&  _  gas: _ 

~  s*  +  IO.ISs  +  G2S 

These  results  compare  favorably  with  information  derived  from  other 
sources.  Calculated  results  derived  for  careful  component  analysis 


6.  3 
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Table  IV 


POSITION  SERVOMECHANISM  CHARACTERISTIC  EQUATION  IDENTIFICATION 
OPEN  LOOP  TRANSFER  FUNCTION  SYSTEM  DYNAMIC  EQUATION 


Model 

Corrected  Model 

Identified  Roots 

10  (40-m  ) 

574 

-4.83  +j  23.45 

c 

s  (s+8+nw) 

s  (s+9.66) 

C  -  .201  0)  =24 

n 

10  (70-m  ) 
c 

s  (s+10-hn,) 
c 

652 

s  (s+10,15) 

-5.075  +j  25 

C  =  .204  (J u  =  25  55 
n 

10  (55-m  ) 
c 

s  (s+10+ra^) 

622 

s  (s+10.24) 

-5.12  +j  24.2 

C  =  .206  <o  =  24.95 
n 

10  (60-m^) 
s  (s+10“hn 

625 

s  (s+10.15) 

-5.075  +j  24.4 

C  =  .208  w  =  25 
n 

76- 


indicated  a  closed  loop  transfer  function  of 


e0  (o 
eL  (s) 


_ S9o _ 

s*  +  (its  t  /.s*)s>  4  590 


6.4 


while  various  open  loop  frequency  response  tests  indicated  the  closed 
loop  response  to  be: 


_  £££122 _ 

s2  +  f/ot  /Js  +  S6s-±  no 

In  view  of  the  non-linearities  involved  in  a  d.c.  servo  motor  it  was 
felt  that  these  model  identified  roots  represented  a  successful 
illustration  of  the  validity  of  the  method. 


Discussion 

These  servomechanism  tests  confirmed  that  the  model  modulation 
method  could  be  successfully  applied  to  obtain  the  governing  operating 
equations  of  a  servomechanism.  Beyond  this  the  tests  emphasized  the 
necessity  of  continuous  excitation  of  the  control  system  for  successful 
model  adaptation.  For  a  linear  control  system  the  model  modulation 
adaptation  will  operate  only  while  the  servo  is  undergoing  excitated 
or  active  response  to  an  input  signal.  It  would  be  expected  that  the 
correction  or  adaptive  signal  would  be  generated  in  response  to  a 
step  for  a  period  equal  to  the  settling  time  of  the  system.  Beyond 
the  settling  time,  with  no  excitation  response  operation  to  generate 
error,  the  time  constants  of  the  adaptive  filter  network  will  allow 
dissipation  of  the  adaptive  signal  magnitude  The  analog  computer 
simulations  detailed  in  this  chapter  illustrate  this  operation  of  the 
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adaptive  loop  For  the  analog  tests  the  general  input  frequency  of 

the  square  wave  was  0.5  cps.  The  average  settling  time  for  the 

second  order  model  was  T  —  4/Su)  ^  10  sec,  and  the  bandwidth  was 

n 

f  ^  1.9  cps  When  tests  were  run  with  systems  with  small  zeta  ,  lower 
input  frequencies  would  achieve  successful  adaptation.  In  all  cases, 
successful  operation  was  attained  with  frequencies  which  were  as  high 
as  the  bandwidth  frequency. 

In  the  case  of  the  operating  servo,  the  non-linearities  of  the 
friction  caused  a  faster  damping  than  the  linear  interpolation  would 
indicate.  For  this  system,  the  adaptation  proceeded  successfully 
when  the  driving  square  wave  frequency  allowed  limited  oscillation 
per  step.  This  in  effect  was  a  partial  masking  of  the  non-linearity 
through  testing  techniques. 

These  situations  emphasize  the  limitation  of  control  system 
excitation  in  the  use  of  the  model  modulation  method  It  is  obvious 
that  this  limitation  exists  for  both  system  adaptation  and  for  system 
identification  through  model  adaptation  A  potential  solution  of 
this  problem  might  be  found  through  investigation  into  utilization 
of  signal  identification  circuits.  This  signal  identification  might 
be  applied  to  the  output  of  the  control  system  to  control  the  operation 
of  the  adaptive  circuit.  This  identifier  could  allow  the  adaptive 
circuit  to  operate  while  the  control  system  would  be  underdoing 
oscillatory  or  excited  response.  When  this  control  system  response 
became  quiescent,  the  identifier  could  place  the  adaptive  loop  into 
a  hold  state  with  a  constant  level  of  output  maintained  to  maintain 
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the  adaptation  until  such  time  as  the  control  system  would  be  again 
activated.  Additionally,  it  is  expected  that  the  adaptive  loop 
elements  would  be  synthesized  through  the  use  of  operational  ampli¬ 
fiers.  Since  these  are  magnitude  sensitive,  another  role  of  the  signal 
identifier  would  be  to  hold  adaptive  loop  operation  when  signal  mag¬ 
nitude  became  excessive. 

The  techniques  of  sampled-data  control  systems  offer  attractive 
possibilities  for  circuits  to  fulfill  these  requirements.  Hybrid 
digital-analog  circuits  might  also  be  investigated.  In  any  case, 
these  techniques  warrant  additional  investigation. 
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Chapter  VIT 


CONCLUSIONS  AND  RECOMMENDATIONS 


Cone lus ions 

It  is  concluded  from  the  results  of  this  investigation  that  the 
identification  of  the  characteristic  equation  roots  of  a  control 
system  by  model  adaptation  through  model  modulation  is  a  valid 
technique.  The  technique  is  subject  to  operating  restrictions  as 
follows : 

1.  The  control  system  must  be  linear  or  be  subjected  to  operating 
non-linearities  of  minor  magnitude. 

2.  The  open  loop  transfer  function  must  have  but  one  or  two 
parameters  whose  magnitude  is  uncertain. 

3.  The  additional  parameters  of  the  open  loop  transfer  function 
must  be  known  or  be  of  sufficient  magnitude  to  be  incon¬ 
sequential. 

4.  The  control  system  must  be  undergoing  active  operation  in 
response  to  an  input  signal* 

It  is  further  concluded  that: 

1.  The  roots  of  the  characteristic  equation  of  a  control  system 
may  be  identified  with  a  high  degree  of  accuracy.  During 
these  tests,  these  dominant  roots  were  identified  to  within 
an  average  error  of  2.2?0  of  the  undamped  natural  frequency 
and  to  within  an  average  error  of  4.76?Q  of  the  complex  phase 
angle . 
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2.  When  the  order  of  the  characteristic  equation  is  unknown, 
second  order  dominant  mode  approximations  will  yield 
consistent  results  which  are  valid  in  the  approximation 

Recommendations 

The  results  of  this  investigation  indicate  that  there  are 
certain  areas  of  interest  which  warrant  further  investigation.  The 
first  of  these  is  investigation  into  the  stability  of  the  dynamic 
adaptive  loop.  The  stability  limits  are  subject  to  various  major 
bounds  in  loop  gain,  input  signal  frequency,  and  input  signal  magnitude. 
Theoretical  analysis  of  these  bounds  will  be  extremely  difficult, 
thus  analog  computer  investigation  might  offer  a  valid  method  of 
investigation. 

The  second  of  these  would  be  the  investigation  of  methods  to 
remove  the  necessity  of  a  continual  excited  response  from  the  control 
system.  Certain  suggestions  have  been  made  in  the  text  which  might 
be  applied  to  achieve  this  end. 

The  third  suggested  area  of  investigation  is  in  the  identifi¬ 
cation  of  control  system  operation  subject  to  known  non-linearities 
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